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Abstract 

It is shown that the chemical equilibrium condition of the system can be unambiguously identified by 
analyzing the isospin evolution of fragments produced in nuclear multifragmentation process. As far 
as the chemical equilibrium is established, the isotope production can be used for finding properties of 
nuclear matter and fragments in the freeze-out volume, for example, via the isoscaling phenomenon. 



1 INTRODUCTION 

Isotopic effects and chemical equilibrium in 
nuclear reactions are receiving increasing at- 
tention because they give a new insight into the 
nuclear matter properties at extreme condi- 
tions. This is of high current interest, in partic- 
ular, for astrophysical applications. Reaction 
of multifragmentation of excited nuclei is very 
promising in this respect since it can be consid- 
ered as manifestation of liquid-gas type phase 
transition in nuclei, therefore, it can directly 
address the nuclear matter behavior at subnu- 
clear densities. A number of theoretical and 
experimental results was already reported on 
this subject [1-6]. Special interest to this prob- 
lem arose after discovering the nuclear caloric 
curve based on the fragment isospin properties 
[J7J and observation of similar isotopic temper- 
atures in many experiments. 

2 STATISTICAL APPROACH 

It is commonly accepted that the detailed 
balance principle is responsible for the chemi- 
cal equilibrium between fragments of different 
sizes. Within statistical models, this can be 
written relationship between their chem- 
ical potentials. In the finite size systems we 
can directly implement this principle by taking 
into account all partitions with corresponding 



weights. For the case of the nuclear multifrag- 
mentation, both methods give very close re- 
sults concerning the fragment isospin. 

The statistical multifragmentation model 
(SMM) is applied for the following analysis. 
This model has been developing since long 
ago and is successful in describing experi- 
mental data (see, e.g., references in || |J). 
Presently, many other statistical models are 
based on similar conceptions, therefore, the 
demonstrated physical results are representa- 
tive for the whole statistical approach. The 
SMM assumes statistical equilibrium of highly 
excited nuclear systems at a low-density freeze- 
out stage p < po/3 (po ~0.15 fm -3 is the normal 
nuclear density) ||. The SMM includes also 
the low energy de-excitation decay modes and 
predicts a natural transition to the multifrag- 
mentation regime from the compound nucleus 
decay ||. All breakup channels (partitions) 
composed of nucleons and excited fragments 
are considered and the conservation of mass, 
charge, momentum and energy is taken into 
account. In the microcanonical treatment the 
statistical weight of decay channel j is given by 
Wj oc exp Sj , where Sj is the entropy of the sys- 
tem in channel j which is a function of the ex- 
citation energy E X) mass number A s , charge Z s 
and other parameters of the source. Different 
breakup partitions are sampled according to 
their statistical weights uniformly in the phase 



space. After breakup, the fragments propagate 
independently in their mutual Coulomb field 
and undergo secondary decays. The deexcita- 
tion of the hot primary fragments proceeds via 
evaporation, fission, or Fermi-breakup [K|. At 
the freeze-out density the SMM treats hot frag- 
ments (A, Z are mass number and charge of 
the fragments) as nuclear matter pieces: Light 
fragments with A < 4 are considered as sta- 
ble particles ("nuclear gas") with masses and 
spins taken from the nuclear tables. Fragments 
with A > 4 are treated as heated nuclear liq- 
uid drops, and their individual free energies 
Faz are parameterized as a sum of the bulk, 
surface, Coulomb and symmetry energy con- 
tributions 

Faz = F^z + f az + e az + e az- (!) 

The standard expressions [§] for these terms 
are: Ff z = (-W - T 2 /e )A, where the 
parameter e is related to the level density, 
and Wo=16 MeV is the binding energy of nu- 

5o^ 2/3 (feS) 5/4 , where 



clear matter; F% z 



B =18 MeV is the surface coefficient, and 
T c =18 MeV is the critical temperature of nu- 
clear matter; E^ z = cZ 2 /A 1 ^ 3 , where c is the 
Coulomb parameter obtained in the Wigner- 
Seitz approximation, c = (3/5)(e 2 /r )(l — 
(/VPo) 1 ^ 3 )) with the charge unit e and r =1.17 
fm; E s A y z = <y{A - 2Z) 2 /A, where 7=25 MeV 
is the symmetry energy parameter. These pa- 
rameters are those of the Bethe-Weizsacker for- 
mula and correspond to the assumption of iso- 
lated fragments with normal density in the 
freeze-out volume, an assumption found to be 
quite successful in many applications. It is 
to be expected, however, that in a more re- 
alistic treatment primary fragments will have 
to be considered not only excited but also ex- 
panded and still subject to a residual nuclear 
interaction between them. These effects can 
be accounted for in the fragment free energies 
by changing the corresponding liquid-drop pa- 
rameters, provided such modifications are also 
indicated by the experimental data. For ex- 
ample, for the symmetry energy, this infor- 
mation may be obtained from the isoscaling 



phenomenon. An unresolved problem for all 
statistical models on the marked is the flow 
development observed at very high excitation 
energies. However, it seems not crucial for the 
isospin studies, since the flow influences mainly 
the kinetic energies of fragments and, possi- 
bly, their charge distributions. As far as the 
charge distributions are described, the isotope 
composition of the fragments can be reliably 
addressed within the statistical approach. 

3 EVOLUTION OF THE 
FRAGMENT ISOSPIN 

The crucial questions to be answered are: 1) 
what is the isospin distribution for fragments 
with fixed A (or Z), 2) how does the isospin of 
fragments change with the isospin of the ther- 
mal source, 3) how does the fragment isospin 
change with excitation energy of the source, 
when the fragment charge distribution evolves. 
In the grand canonical approximation ||, the 
mean multiplicity of fragments is given by 



(N AZ ) 



gAz-r^Ai exp 



Faz - \lA — vZ 



T 



(2) 

where qaz is the degeneracy factor, At is the 
nucleon thermal wavelength, Vf is the "free" 
volume, and fi and v are the chemical poten- 
tials responsible for the mass and charge con- 
servation in the system, respectively || [|. As 
shown in Ref. [i~0|, the average charge (Z a) of 
fragments with mass A and the width a z of 
their charge distribution can be written as 



(Za) 



(47 + v)A 
87 + 2cA 2 / 3 ' 




The symmetry energy, Coulomb interaction 
and the chemical potential v are directly re- 
sponsible for the isospin of the produced frag- 
ments. As seen from Eq. (||), the neutron con- 
tent of the fragments increases with their mass 
number, because of the Coulomb term. This 
trend exists in the stable nuclei and it presents 
also in multifragmentation experimental data. 
One can see it from Fig. 1, where the neutron 
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Figure 1: Measured N/Z (top) and mean value 
(N/Z) (bottom) as a function of the charge Z of 
projectile fragments from the reaction 197 Au + Cu 
at E/A=600 MeV. The lines in the top panel 
represent the valley of stability according to the 
Weizsacker mass formula (long dashed) and EPAX 



parameterization [11] (short dashed). In the bot- 
tom panel, the experimental data are given by the 
open circles and the prediction of the SMM by the 
histogram. 

to proton (N/Z) ratio of fragments observed by 
ALADIN |L2] at multifragmentation of Au-like 
sources is shown. The agreement of the SMM 
calculations with the data is obtained simul- 
taneously with a very good description of the 
fragment partitions [|T^], this is an important 
requirement for this kind of analysis. 

An intuitive expectation that the fragments 
produced from a neutron rich source are also 
neutron rich was confirmed in many experi- 
ments (see, e.g., [||, |5J). It is fully consis- 
tent with the statistical predictions: As cal- 
culations show ||, |5| the chemical potential 
v decreases with the neutron excess of the 
source, therefore, this effect for intermediate 
mass fragments (IMF, charges Z=3-20) is ex- 
plained by Eq. @. 
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Figure 2: The neutron-to-proton ratio N/Z (top) 
and relative yield (bottom) of hot primary frag- 
ments produced after the break-up of Au nuclei 
at different excitation energies: 3 (solid lines), 4 
(dashed lines), 5 (dotted lines) and 8 (dot-dashed 
lines) MeV per nucleon. 



In Fig. 2 the N/Z ratio of the fragments to- 
gether with their mass distributions are shown. 
The calculations were done with the micro- 
canonical Markov-chain SMM version || to 
take into account the mass and charge con- 
servation explicitly. The results are consistent 
with the grand canonical ones (see also || |J), 
however, one can see that the finite-size effects 
favor the saturation of the N/Z ratio of very 
large fragments with A > A s /2. This deviation 
from Eq. (|3]) is important for low multiplicity 
channels at small E x . 

The SMM predicts an interesting behavior 
of fragment isospin with changing fragment 
charge distribution. In Fig. 2 one can also 
see the evolution of the N/Z ratio and mass 
distribution of fragments in the excitation en- 
ergy range E x =3-8 MeV/nucleon. It is seen 
that the fragment mass distribution evolves 
from the U-shape, at the multifragmentation 
threshold E x =E th 3 MeV/nucleon, to an ex- 
ponential fall-off at high energies. During this 



evolution the temperature reaches a "plateau" 
and is nearly constant J7], |[]. As the energy- 
increases the N/Z ratio of primary IMFs' in- 
creases, too. The reason is that the heavi- 
est neutron-rich fragments are destroyed at in- 
creasing excitation energy, and most of their 
neutrons are bound in IMFs, since the number 
of free neutrons is still small at this stage. Gen- 
erally, the N/Z ratios of fragments decrease af- 
ter their secondary deexcitation (e.g., compare 
Fig. 1 and Fig. 2). If the temperature increases 
considerably with the excitation energy (in the 
energy region of small E x < E t h) the N/Z ra- 
tios of final cold fragments can decrease with 
the energy. However, when the temperature 
reaches the "plateau" the internal excitation 
energy of hot fragments does not change (see 
Eq. (HD) and the final isotope composition be- 
comes proportional to the primary fragment 
isospin. In this case, the effect of increasing the 
N/Z ratio of IMFs can survive after the sec- 
ondary deexcitation. The experimental data 
H show that this effect exists and is consistent 
with the SMM predictions. Therefore, it is a 
good tool for probing the freeze-out conditions 
and verification of the chemical equilibrium. 

It is important that the demonstrated isospin 
evolution takes place in the energy range which 
is usually associated with a liquid-gas type 
phase transition in finite nuclei ||. The re- 
alistic statistical model, designed to reproduce 
the experiment, does not support conclusions 
of some mean- field approaches (see Ref. fl3l ) 
about separation of the nuclear matter into a 
neutron-rich gas and isospin symmetric "liq- 
uid" matter during the phase transition. On 
the contrary, the chemical equilibrium in fi- 
nite nuclear systems suggests few free neutrons 
and neutron-rich big fragments at the freeze- 
out stage. That is consistent with experimen- 
tal observations. 

4 ISOSCALING 

By using the chemical equilibrium conception 
guideline one can extract more informa- 
tion about fragments and nuclear matter prop- 



erties from the analysis of experimental data. 
The scaling properties of cross sections for frag- 
ment production with respect to the isotopic 
composition of the emitting systems were stud- 
ied long ago in light-ion induced reactions 0. 
Recently, the isoscaling was also reported for 
heavy- ion reactions |j| . It is constituted by the 
exponential dependence of the production ra- 
tios R21 for fragments with neutron number N 
and charge Z in reactions with different isospin 
asymmetry: 



R 



21 



Y 2 (N, Z) 
Y 1 (N,Z) 



C-exp(N -a + Z-P), (4) 



with three parameters C, a and (3. Here Y 2 and 
Yi denote the yields from the more neutron rich 
and the more neutron poor reaction system, 
respectively. In some reactions, the parameters 
a and /3 have the tendency to be quite similar 
in absolute magnitude but of opposite sign 0. 
This suggests an approximate scaling with the 
third component of the isospin t 3 =(N — Z)/2: 



R 



12 



Yi(N,Z) 
Y 2 {N, Z) 



C ■ exp(-t 3 ■ # 



(5) 



As an example, the isotope ratios measured in 
five pairs of reactions induced by light particles 
from protons of 0.66 GeV to a-particles of 15.3 
GeV incident energy on 112 Sn and 124 Sn tar- 
gets, and normalized with respect to the ratio 
for 6 Li, are shown in Fig. 3 (see Ref. ||). The 
parameter f3 t3 decreases from 1.08 to 0.68 with 
increasing energy. The inverse scaling parame- 
ter 1 / p t 3 is found to increase approximately in 
proportion to the isotope temperature deduced 
from the yields of He and Li isotopes. 

The isoscaling phenomenon arises naturally 
in a statistical fragmentation mechanism |], 
14]. As evident from Eq. (0), for two sys- 
tems 1 and 2 with different total mass num- 
bers (Ai and A 2 ) and charges {Z\ and Z 2 ) but 
with the same temperature and density, the 
ratio of fragment yields produced in these sys- 
tems is given by Eq. @ with parameters a = 
(//1 —fJv)/T and (3 = ((/xj. -/z 2 ) + [y x -v 2 ))/T. 
In the SMM the isoscaling parameters reflect 
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Figure 3: Normalized isotopic effect R±2 for five 
reactions, labeled with the projectile energy (in 
GeV) and shifted by factors of three for clarity. 
The symbols are for H, He, Li, Be, and B frag- 
ments. The lines are the results of the exponential 
fits according to Eq. (||). 

properties of the produced fragments, in par- 
ticular, the widths of their charge (mass) distri- 
butions (see Eq. (|3])). As shown in Ref. || the 
potential differences depend essentially only on 
the coefficient 7 of the symmetry term and on 
the isotopic compositions of the sources but 
not on the temperature: 



/xi-/x 2 « -47(^2 -342). 



(6) 



This formula is a very good approximation 
(within 3%) in the grand-canonics, and it was 
confirmed for the multifragmentation region 
(T>5 MeV) with the microcanonical Markov 
chain SMM calculations. Moreover, since the 
temperature can be found in an independent 
way, e.g. as the isotope temperature, Eq. @ 
provides an effective way for determination of 
the symmetry energy coefficient 7 via experi- 
mental a and (3 parameters. 




Figure 4: Ratios of isotope yields produced at the 
breakup of 112 Sn and 124 Sn sources from Markov- 
chain SMM calculations for three excitation ener- 
gies E x /A = 2, 5 and 8 MeV. The top and bottom 
panels are for hot and cold fragments, respectively. 
The solid lines give the logarithmic slope param- 
eters /3j3=1.74, 1.23, 1.09 for hot fragments, and 
1.33, 0.95, 0.84 for cold fragments. 



A typical isoscaling obtained with the SMM 
is shown in Fig. 4. These calculations can be 
used to estimate quantitatively the influence of 
the secondary deexcitation on the observable 
parameters. The isoscaling parameters can 
change, and the corresponding correction to 
the extracted symmetry energy coefficient has 
to be applied. The full analysis was performed 
in Ref. |J and the value 7 ~ 22.5 MeV was 
reported, which has a tendency to be smaller 
than the conventional value. Provided it can 
be substantiated by other data and analyses, 
this would indicate that the symmetry part of 
the fragment binding energy is slightly weaker 
than that of isolated nuclei. Fragments, as 
they are formed at breakup, may have a lower 
than normal density. 



5 CONCLUSION 

Conception of chemical equilibrium is very 
fruitful for the theoretical description of multi- 
fragmentation. The chemical equilibrium can 
be unambiguously verified by analyzing isospin 
content of the fragments and its evolution to- 
gether with the fragment yields. A convincing 
proof of this condition is an excellent descrip- 
tion of the fragment charge partitions and the 
isotope distribution of fragments obtained in 
different experiments within statistical models 
A remarkable trend char- 



12 



(e.g., see 

acterizing the chemical equilibrium is an in- 
crease of the N/Z ratio with the mass num- 
ber of fragments f|, [10| . This is quite instruc- 
tive because dynamical processes can lead to 
the opposite ("decreasing") trend, as shown 
by M. Colonna and M.Di Toro in 0. Ex- 
perimental observations of this trend may be 
obscured by the secondary deexcitation. How- 
ever, by studying evolution of the N/Z ratio 
of fragments with changing their charge dis- 
tribution one can identify it for sure: The ef- 
fect of increasing N/ Z of IMFs within the en- 
ergy range of the phase transition, observed 
for the Au sources ||, can be explained in the 
case of the " increasing" trend only. The chem- 
ical equilibrium conception can also provide a 
key for explanation of many other phenomena. 
For example, as shown in |4j], the midrapid- 
ity emission of neutron-rich IMF in peripheral 
heavy-ion collisions can be explained as a re- 
sult of the Coulomb proximity of the sources 
and the angular momentum influence within 
the statistical picture of the two (projectile- 
like and target-like) decaying sources. An im- 
portant advantage of this conception is that it 
allows for obtaining direct experimental infor- 
mation about nuclear properties at the break- 
up. A well-known example is the isotope tem- 
peratures, extracted from experimental data 
0. Another example is the isoscaling phe- 
nomenon, which is connected with the symme- 
try energy of nuclear matter. As shown here, 
the isoscaling gives a new way for probing the 
symmetry energy term of the fragments in the 



freeze-out volume. All these observations can 
be directly used in the astrophysical cases (su- 
pernovas, neutron stars), where the chemical 
equilibrium, is believed, to be established. 
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this presentation. 
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